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Abstract 

A  computational  fuel  cell  dynamics  framework  is  used  to  develop  a  unified  water  transport  equation  for  a  proton  exchange  membrane  fuel 
cell  (PEMFC).  Various  modes  of  water  transport,  i.e.,  diffusion,  convection  and  electro-osmotic  drag,  are  incorporated  in  the  unified  water 
transport  equation.  The  water  transport  model  is  then  applied  to  elucidate  water  management  in  three-dimensional  fuel  cells  with  dry-to-low 
humidified  inlet  gases  after  its  validation  against  available  experimental  data  for  dry  oxidant  and  fuel  streams.  An  internal  circulation  of  water 
with  the  aid  of  counter-flow  design  is  found  to  be  of  vital  importance  for  low-humidity  operation,  for  example,  in  the  portable  application  of 
a  PEMFC  without  an  external  humidifier.  The  general  features  of  water  transport  in  PEMFCs  are  discussed  to  show  various  water  transport 
regimes  of  practical  interest,  such  as  anode  water  loss,  cathode  flooding,  and  the  equilibrium  condition  of  water  at  the  channel  outlets, 
particularly  for  limiting  situations  where  anode  and  cathode  water  profiles  acquire  an  equilibrium  state.  From  the  practical  point  of  view, 
the  effects  of  the  flow  arrangement,  membrane  thickness,  and  inlet  gas  humidity  as  important  determinants  of  fuel  cell  performance  are  also 
analyzed  to  elucidate  fuel  cell  water  transport  characteristics. 

©  2005  Elsevier  B.V.  All  rights  reserved. 

Keywords:  Proton  exchange  membrane  fuel  cell;  Computational  study;  Water  transport 


1.  Introduction 

In  the  previous  modelling  efforts  of  proton  exchange  mem¬ 
brane  fuel  cells  (PEMFCs)  [1-4],  the  ionomer  membrane 
remained  fully  humidified,  which  yields  a  constant  proton 
conductivity  and  minimum  ohmic  loss.  Under  such  oper¬ 
ating  conditions,  however,  the  gas  phase  is  over- saturated 
with  water  vapour  and  water  condensation  may  occur  on 
the  cathode  side  even  at  low  operating  current  densities.  The 
ensuing  cathode  flooding  by  liquid  water  subsequently  hin¬ 
ders  the  access  of  oxygen  to  the  cathode  catalyst  layer  and 
results  in  significant  concentration  polarization.  In  addition, 
full  humidification  of  the  reactant  gases  externally  involves 
parasitic  volume  and  power  losses  in  the  fuel  cell  system. 
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Therefore,  it  is  of  practical  interest  to  operate  PEMFCs 
under  low-humidity  conditions,  while  still  maintaining  the 
polymer  electrolyte  membrane  adequately  hydrated.  Such 
low-humidity  operating  strategies  will  also  reduce  cathode 
flooding  and  space  requirement  for  auxiliary  equipment  to 
a  minimum  extent.  Obviously,  such  an  innovative  operating 
scheme  hinges  upon  a  good  understanding  of  water  genera¬ 
tion,  transport,  and  distribution  within  PEMFCs. 

The  research  presented  here  is  dedicated  to  studying  water 
transport  and  distribution  in  PEMFCs  by  extending  the  com¬ 
putational  fuel  cell  dynamics  (CFCD)  model  developed  by 
Um  et  al.  [5].  A  key  challenge  in  applying  a  unified  CFCD 
framework  to  describe  water  transport  is  to  recognize  the 
fact  that  there  are  different  phases  of  water  existing  in  var¬ 
ious  regions  of  a  PEMFC,  i.e.  water  in  the  gas  phase  and 
in  the  membrane.  As  a  result,  phase  equilibrium  must  be 
considered.  Furthermore,  various  modes  of  water  transport, 
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Nomenclature 

a  effective  catalyst  area  per  unit  volume 

(cm2  cm-3)  or  water  activity 
c  molar  concentration  (mol  cm-3) 

D  mass  diffusivity  of  species  (cm2  s-1) 

EW  equivalent  weight  of  membrane 

I  current  density  (A  cm-2) 

j  transfer  current  (A  cm-3) 

rid  electro-osmotic  drag  coefficient 

P  pressure  (Pa) 

R  gas  constant,  8.3 14  J  mol-1  K-1 

t  time  (s) 

T  temperature  (K) 

u  velocity  vector  (cm  s  _  1 ) 

U  inlet  velocity  at  gas  channel 

V  potential  (V) 

y  axial  location 

Greek  letters 

f  stoichiometric  flow  ratio  (Eq.  (1)) 

k  ionic  conductivity  (S  m_  1 ) 

p  density  (kg  cm-3) 

Subscripts 
a  anode 

c  cathode 

crit  critical  water  loss  in  the  anode 

e  electrolyte 

fl  cathode  flooding 

g  gas  diffusion  layer 

ref  reference  value 

W  water 

Superscripts 
eff  effective  value 

g  gas  phase 

m  membrane  phase 

sat  saturation  value 


i.e.,  diffusion,  convection  and  electro-osmotic  drag,  must  be 
incorporated  in  the  unified  water  transport  equation.  An  equi¬ 
librium  water  uptake  curve  in  the  membrane  phase  [1]  is  thus 
introduced  to  extend  the  existing  CFCD  framework  to  various 
low  membrane  humidification  conditions. 

Water  is  produced  at  the  cathode  catalyst  layer  as  a  result  of 
the  oxygen  reduction  reaction  (the  product  water)  and  can  be 
brought  into  the  cell  by  humidified  reactant  gases.  The  various 
modes  of  water  transport  through  the  membrane  electrolyte, 
primarily  the  electro-osmotic  drag  effect  and  diffusion  driven 
by  the  concentration  gradient,  are  illustrated  in  Fig.  E  The 
water  flux  due  to  electro-osmotic  drag  is  directly  propor¬ 
tional  to  the  proton  flux  (/ceiiAF);  the  proportionality  constant 
is  called  the  ‘electro-osmotic  drag  coefficient’,  rid.  The  dif¬ 
fusion  flux  of  water  in  the  membrane  is  usually  described 


by  a  water  diffusion  coefficient  and  by  the  gradient  of  molar 
concentration  of  water. 

A  number  of  excellent  computational  efforts  on  water 
management  in  PEMFCs  have  been  made  [6-10]  to  account 
for  multiple  water  transport  modes  across  die  proton 
exchange  membrane.  Some  approaches  solve  a  simplified 
one-dimensional  water  transport  model  [6,7,11],  other  com¬ 
putational  works  excluded  the  thin  catalyst  layer  or  treat  it  as 
an  interface  [8-10, 12] .  A  significant  contribution  to  analysing 
water  transport  in  PEMFCs  has  been  made  by  Zawodzinski 
et  al.  [13-15]  and  Springer  et  al.  [1,16],  both  experimentally 
and  computationally. 

Zawodzinski  et  al.  [13-15]  conducted  a  series  of  experi¬ 
ments  on  the  sorption  characteristics  of  Nation®  membrane. 
They  measured  water  uptake  by  a  membrane  in  equilibrium 
with  liquid  water  and  water  vapour,  and  estimated  die  electro- 
osmotic  drag  coefficient  for  various  membrane  water  con¬ 
tents.  A  coefficient  of  unity  for  membranes  equilibrated  with 
water  vapour  over  a  wide  range  of  water  activities  was  deter¬ 
mined  experimentally.  Ren  and  Gottesfeld  [17]  showed  the 
effects  of  temperature,  membrane  hydration  level,  membrane 
equivalent  weight  and  membrane  type  on  the  electro-osmotic 
drag  coefficient.  They  found  that  the  drag  coefficient  is  inde¬ 
pendent  of  the  cell  current  density. 

A  diffusion  model  was  introduced  by  Springer  et  al.  [1] 
for  a  Nation®  117  membrane,  which  included  the  electro- 
osmotic  drag  term.  They  assumed  an  equilibrium  between  the 
gas  phase  and  the  membrane  phase  of  water  in  the  Nation® 
membrane,  and  then  determined  the  membrane  water  con¬ 
tent  at  the  interface  from  the  activity  of  water  vapour.  Hsing 
and  Futerko  [9]  developed  a  two-dimensional,  finite  element 
based  model  for  a  PEMFC  without  external  humidifica¬ 
tion  of  the  inlet  gas  streams.  The  flow  field  was  approxi¬ 
mated  using  a  potential  flow  equation  with  a  stream  func¬ 
tion.  The  membrane  water  content  at  the  anode  vertical  line 
membrane  interface  was  calculated  under  the  assumption  of 
phase  equilibrium  between  the  vapour  and  the  membrane 
phases,  and  the  cathode  vertical  line  membrane  interface 
was  taken  to  be  fully  hydrated.  More  recently,  Janssen  [10] 
demonstrated  one-  and  two-dimensional  numerical  models  to 
account  for  water  transport  through  the  membrane  electrode 
assembly  (MEA)  under  various  humidification  conditions  of 
inlet  gases.  Computed  results  were  compared  with  exper¬ 
imental  data  measured  by  Janssen  and  Overvelde  [18].  In 
this  model,  it  is  assumed  that  the  cell  current  density  is  con¬ 
stant  everywhere,  and  the  catalyst  layers  are  excluded  in  the 
water  transport  calculation.  Reacting  species  along  the  chan¬ 
nel  and  in  the  through-membrane  direction  of  the  MEA  are 
calculated  based  on  the  assumption  of  one-dimensional  mass 
transport. 

In  the  present  work,  a  unified  water  transport  equation 
within  the  single-domain  CFCD  framework  is  developed 
with  the  aid  of  the  equilibrium  water  uptake  curve  in  the 
membrane  phase.  Emphasis  is  placed  on  how  to  derive 
correctly  the  unified  water  transport  equation  with  the 
various  transport  modes  incorporated.  Also,  this  water 
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Fig.  1.  Schematic  diagram  of  water  transport  mechanism  in  fuel  cell  domain. 


transport  model  is  validated  against  available  experimental 
data  from  the  literature,  and  then  applied  to  elucidate 
water  management  in  three-dimensional  fuel  cells  with  low 
humidified  inlet  reactant  gases. 

2.  Characteristics  of  water  concentration  profiles  in 
anode  and  cathode 

The  general  features  of  water  transport  in  a  fuel  cell 
domain  are  schematically  shown  in  Fig.  1.  Water  leaves 
the  anode  under  the  electro- osmotic  drag  when  an  electri¬ 
cal  potential  is  applied,  while  water  is  driven  back  from 
the  cathode  to  the  anode  by  the  concentration  gradient  that 
primarily  arises  from  water  production  at  the  cathode  cata¬ 
lyst  layer.  Water  loss  on  the  anode  surface  of  the  membrane 
caused  by  the  electro-osmotic  drag  becomes  significant  at 
large  current  densities  and  the  back-diffusion  of  water  usu¬ 
ally  becomes  stronger  along  the  axial  direction  as  depicted  in 
Fig.  1. 

Typical  water  concentration  profiles  in  the  anode  and  cath¬ 
ode  along  the  flow  direction  are  depicted  in  Fig.  2.  The  water 
concentration  at  the  anode  vertical  line  membrane  interface 
decreases  along  the  axial  direction  and  then  reaches  the  min¬ 
imal  point  at  which  the  net  water  flux  across  the  membrane 
becomes  zero.  This  implies  that  the  electro-osmotic  drag  is 
balanced  by  the  back-diffusion  of  water.  The  location,  ycrit 
in  Fig.  2  indicates  this  zero  net  water  flux  location.  Before 
this  critical  location,  the  anode  is  characterized  by  water  loss, 
while  after  the  location  the  anode  gains  water.  In  the  worst  sit¬ 
uation,  the  anode  can  be  fully  dehydrated  if  the  back-diffusion 


from  the  cathode  is  not  sufficiently  strong  to  compensate  for 
the  loss  of  anode  water  by  electro-osmotic  drag.  This  critical 
location  is  also  important  in  that  it  signifies  the  most  resistive 
location  of  the  membrane. 

For  current  densities  greater  than  a  certain  value,  so  much 
water  is  produced  at  the  cathode  that  the  gas  is  saturated 
with  water  vapour  and  a  vapour-to-liquid  phase  change  is 
imminent,  which  brings  about  cathode  flooding.  Theoreti¬ 
cally,  liquid  water  appears  at  a  flooding  point,  ya,  shown  in 
Fig.  2. 

When  a  thin  membrane  is  used  in  a  fuel  cell,  the  water 
exchange  between  the  anode  and  cathode  is  greatly  facili¬ 
tated.  Thereby,  the  water  concentrations  on  the  anode  and 
cathode  sides  of  the  membrane  approach  each  other  towards 
the  exit,  a  state  called  the  ‘equilibrium  of  water’.  This  is 
depicted  as  curve  Cl  in  Fig.  2.  This  equilibrium  condition 
occurs  only  when  the  back-diffusion  of  water  through  the 
membrane  is  largely  dominating.  If,  however,  the  membrane 
is  sufficiently  thick  that  back-diffusion  of  water  cannot  com¬ 
pensate  for  the  electro-osmotic  drag  of  water  from  the  anode 
to  cathode,  the  differential  in  water  concentration  between 
the  anode  and  the  cathode  will  widen  as  a  result  of  water 
production  on  the  cathode  side.  This  is  shown  as  curve  C2  in 
Fig.  2. 

3.  CFCD  model  for  water  transport 

For  this  water  transport  study,  a  multi-dimensional  CFCD 
model  developed  by  Um  et  al.  [5]  is  extended  to  model  accu¬ 
rately  water  transport  and  distribution  in  a  PEMFC,  as  shown 
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Fig.  2.  Graphical  representation  of  water  concentration  profiles  along  flow  direction  in  both  anode  and  cathode  of  co-flow  design,  showing  various  regimes  of 
water  transport. 


in  Fig.  3,  particularly  for  low-humidity  operation  in  which  the 
anode  and  the  cathode  inlet  gases  are  only  slightly  humidi¬ 
fied. 

A  distinguishing  feature  of  the  present  CFCD  model  lies  in 
its  unified  treatment  of  water  transport  throughout  all  regions 
of  a  fuel  cell.  The  model  how  the  single  species  equation 
for  water  presented  in  the  previous  work  [5]  can  describe  all 
transport  modes  in  various  regions. 

3.1.  Flow  channels  and  backing  layers 

When  applied  to  water,  the  species  conservation  equation 
for  water  reduces  to  the  following  for  flow  channels  and  back¬ 


ing  layers: 

+  V  .  (euc%)  =  V  •  (D$fv4)  (1) 

where  s  is  the  porosity  of  the  gas  channel  and  the  backing 
layer;  c ^  the  water  molar  concentration  in  gas  phase;  u  the 
intrinsic  fluid  velocity  and  D ^  is  an  effective  mass  diffusivity 
defined  in  the  previous  study  [5]. 

The  physical  meaning  of  Eq.  (1)  is  apparent,  each  term 
has  a  unit  of  (molm-3  s-1).  An  effective  transport  property 
for  the  porous  region  is  modified  from  its  intrinsic  physical 
property  by  a  Bruggemann  factor  [5]. 


Fig.  3.  Description  of  cell  geometry  of  PEMFC. 
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3.2.  Membrane 

The  general  water  transport  equation  reduces  to  the  fol¬ 
lowing: 


in  order  to  obtain  the  physically  meaningful  variable  more 
conveniently. 

The  diffusion  coefficient  for  in  the  membrane  region 
is  thus  equal  to: 


=  V  •  (Z4v4)  -  f  V  •  (nd/e)  (2) 

for  the  membrane  region.  Here,  use  has  been  made  of  the 
assumption  that  fluid  velocity  is  zero  in  the  membrane.  In 
Eq.  (2),  the  water  molar  concentration,  c^,  is  a  fictitious 
vapour  concentration  (c|^  =  (Psat//?rceii)tf)  that  is  in  ther¬ 
modynamic  equilibrium  with  a  membrane  of  water  content 
X  (or  its  equivalent  membrane  water  concentration,  c ^  = 
(/Odry/EW)A).  The  one-to-one  correspondence  between 
and  c™  can  be  called  the  ‘equilibrium  uptake  curve’ ,  as  shown 
in  Fig.  4(a). 

The  main  advantage  of  using  the  fictitious  gas-phase  water 
concentration,  c^,  in  the  membrane  region  is  to  ensure  con¬ 
tinuity  of  the  variable  field  across  the  electrode  vertical  line 
membrane  interface,  as  shown  in  Fig.  4(b).  The  solid  line  rep¬ 
resents  a  physical  variable,  whereas  the  dashed  line  denotes 
a  fictitious  quantity  that  is  solved  in  the  numerical  program 


(a) 


(b) 

Fig.  4.  (a)  Equilibrium  uptake  curves  for  Nation®  membrane  relating  mem¬ 
brane  water  concentration  to  equilibrium  gas  phase  water  concentration 
(created  from  Springer  et  al.  [1]  relation  between  the  water  content  and  activ¬ 
ity)  and  (b)  schematic  profiles  of  water  concentrations  in  various  regions  of 
a  PEMFC. 


D^  =  D 


c\rm 
m  U6W 

wvy 
d  c 


w 


1  RT 

EW^ 


dX 

da 


where  the  derivatives,  dp/dX  and  dX/da ,  can  be  evaluated 
from  p-X  and  X-a  property  correlations.  Substituting  the 
expression  for  Eq.  (3),  into  Eq.  (2)  results  in  Springer’s 
model  [1]  of  water  transport  through  the  membrane  under 
steady-state: 

V-(D^Vc^)-4v-(nd/e)  =  0  (4) 

r 

or 


V-  (z^^Va)  -  2v.(nd/e)  =  0  (5) 

if  the  membrane  density  is  assumed  to  be  constant  (i.e.,  no 
membrane  swelling  is  considered). 

Therefore,  the  water  transport  equation  used  in  the  present 
model  unifies  the  water  transport  process  in  the  gas  phase  with 
that  in  the  membrane,  and  thereby  permits  a  single-domain 
simplification  in  computations.  A  salient  feature  of  this  uni¬ 
fied  water  transport  equation  is  the  fact  that  the  water  flux 
intrinsically  remains  continuous  across  the  interface  between 
the  gas-diffusion  electrode  and  the  membrane.  This  water  flux 
continuity  is  not  automatically  guaranteed  in  many  previous 
models  using  a  multi-domain  approach. 

A  comment  is  in  order  on  the  effect  of  the  electro- osmotic 
drag  on  water  transport  through  the  membrane  phase.  It  fol¬ 
lows  from  Eq.  (5)  that  this  effect  is  directly  proportional  to 
the  gradient  of  the  ionic  current,  V*/e. 

Various  expressions  of  the  electro-osmotic  drag  coeffi¬ 
cient,  nd,  have  been  proposed  in  the  literature;  for  example, 

'  2.5X 
~22~ 

«d  =  <  1  ,  1 

U0.35A)4  1.474 

1 
V 

The  electro-osmotic  drag  coefficient  expression  given  by 
Springer  et  al.  [1]  has  a  linear  dependence  on  the  water  con¬ 
tent,  as  expressed  by  Eq.  (6a),  whereas  Fuller  [19]  determined 
the  water  drag  coefficient  by  Eq.  (6b).  On  the  other  hand,  the 
experimentally  measured  proton  transport  number  of  water 
shows  a  constant  value  close  to  1.0  for  a  water  content  level 
up  to  14,  while  for  a  membrane  immersed  in  liquid  water, 
the  water  drag  by  protons  has  an  abrupt  discontinuity  and 
jumps  to  2.5  H20/H+  [15].  In  the  present  study,  the  constant 
electro-osmotic  drag  coefficient  given  in  Eq.  (6c)  is  used  for 
computational  purpose. 
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Table  1 


Water  transport  equation  in  various  regions  of  fuel  cell 


Components 

Governing  equation 

Diffusivity 

Anode  gas  channel 

Dg 

Anode  backing  layer 

L5  Dg 

Anode  catalyst  layer 

Slc5  Z4  +  (6mf 

Membrane 

^  +  V  •  (£MCw) 

Dm 

usn 

=  V  •  (DeffVcw)  +  Sw 

Cathode  catalyst  layer 

e[-5  D%  +  D™, 

Cathode  backing  layer 

L5  Dg 
fcg 

Cathode  gas  channel 

Dg 

of  the  numerical  solution  procedure  and  the  code  have  been 
given  in  previous  work  [5].  Stringent  numerical  tests  were 
performed  to  ensure  that  the  solutions  were  independent  of 
the  grid  size.  A  50  x  100  x  60  mesh  (or,  about  300,000  ele¬ 
ments)  was  found  to  provide  sufficient  spatial  resolution.  The 
coupled  set  of  equations  was  solved  iteratively,  and  the  solu¬ 
tion  was  considered  to  be  convergent  when  the  relative  error 
in  each  field  between  two  consecutive  iterations  was  less  than 

icr5. 


4.  Results  and  discussion 


3.3.  Catalyst  layers 


Water  transport  is  through  both  gas  and  membrane  phases 
in  catalyst  layers.  Assuming  that  the  water  diffusion  paths 
through  the  gas  and  membrane  phases  are  in  parallel,  the 
total  diffusive  flux  as  follows: 

drm 

Vw  =  £*-5(-£>wV<4)  -  (£m)"£>WTirV4  =  -£>w  Vt& 

dcw 

(7) 


where 

ne ff  _  U-5  n&  _i_  n 
lJW  ~  £c  lJW  T-  Lf- 


m  U6W 

WTX 
d  c 


w 


Note  that  water  diffusion  is  practically  through  the  gas 
phase  in  the  catalyst  layer  because  the  water  diffusion  coef¬ 
ficient  in  the  membrane  phase  is  roughly  two  orders  of  mag¬ 
nitude  smaller  than  the  gas  phase. 

The  unified  water  transport  equation  is  summarized  in 
Table  1  along  with  its  diffusivity  expressions  in  various 
regions  of  the  PEMFC.  More  information  on  general  govern¬ 
ing  equations  and  source/sink  terms  for  fuel  cell  modelling 
is  given  by  Um  et  al.  [5]. 


This  section  focuses  on  applying  the  unified  water  trans¬ 
port  model  to  three-dimensional  fuel  cells.  This  enables  the 
investigation  of  water  transport  phenomena  and  its  effect  on 
the  cell  performance.  A  series  of  three-dimensional  simula¬ 
tions  are  performed,  including  both  anode  and  cathode  gas 
streams  with  full  humidification  (AFCF),  fully  humidified 
anode  flow  but  dry  cathode  flow  (AFCD),  both  anode  and 
cathode  with  low  humidification  (i.e.,  RH  =  20%)  (AFCF), 
and  both  anode  and  cathode  with  dry  gases  (ADCD).  The 
focus  of  these  simulations  is  on  the  prediction  of  water  con¬ 
tent  profiles  across  the  ionomer  membrane  as  well  as  the 
resulting  current  density  distribution. 

The  feed  stream  in  the  anode  gas  channel  is  assumed 
to  consist  of  hydrogen  and  water  vapour  only  to  simulate 
direct-hydrogen  fuel  cell  vehicles,  while  air  is  supplied  to  the 
cathode  (i.e.,  oxygen,  nitrogen,  and  water  vapour). 

It  is  assumed  that  the  amount  of  Pt  catalyst  and  its  trans¬ 
fer  current  density  within  the  catalyst  layer  remain  the  same 
in  all  simulation  cases.  The  physical  and  transport  prop¬ 
erties  and  the  electrochemical  kinetics  are  summarized  in 
Tables  2  and  3.  The  basic  operating  conditions  and  the  cell 
design  parameters  are  summarized  in  Tables  4  and  5. 


3.4.  Boundary  conditions 


4.1.  Experimental  validation 


At  the  fuel  and  oxidant  inlets,  the  following  conditions  are 
prescribed: 

^  in, anode  =  Uajn,  win,  cathode  =  f^c,in? 

CH2,  anode  —  CH2,a>  CO2,  cathode  —  C02,c? 

anode  —  cH20,a?  CH20, cathode  —  cH20,e  (9) 

The  inlet  velocities  of  fuel  and  oxidant  can  also  be  expressed 
by  their  respective  stoichiometric  flow  ratios,  i.e.,  £a  and  at 
1  A  cm-2 .  At  the  outlets,  both  anode  and  cathode  channels  are 
assumed  to  be  sufficiently  long  so  that  velocity  and  species 
concentration  fields  are  fully  developed. 

3.5.  Numerical  procedures 

The  conservation  equations  [5]  were  discretized  using  a 
finite  volume  method  [20]  and  were  solved  using  a  general- 
purpose  computational  fluid  dynamic  (CFD)  code.  Details 


To  validate  the  unified  water  transport  model  developed 
here,  a  special  application  of  fuel  cells  operated  under  no 
humidity  is  considered.  The  computed  polarization  curve  and 
the  water  removal  rate  through  the  anode  channel  are  com¬ 
pared  in  Fig.  5  with  the  published  experimental  data  of  Biichi 


Table  2 

Electrochemical  kinetics 


Description 

Anode 

Cathode 

Reference  cell  potential, 

0 

Ucii 

<Z>s  (V) 

Equilibrium  potential,  Vo 

o 

o 

II 

O 

Vc,oc  =  0.00257+0.2329 

(V)  [22] 

Exchange  current 

1.0  X  109 

2  x  104  exp[0.014189(T  —  353)] 

density  x  reaction 
surface  area,  a 
(A  m-3) 

Transfer  coefficient,  a 

aa  +  ac  =  2 

ac  =  1 

Faraday  constant,  F  (C) 

96487 

Table  3 

Transport  properties 
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Property 

Value 

Reference 

Ionic  conductivity  of  membrane,  k  (S  m“ 1 ) 

(0.5139A  —  0.326)  exp 

[l268  (  303  7)] 

[1] 

H2O  diffusion  coefficient  in  membrane,  D ^  (m2  s-1) 

/  3.1  X  10“7A(e°-28A  -  1  )e-2346/r 
w  _  1  4.17  x  1CT8A(1  +  161e"7)e“2346/r 

forO  <  A.  <  3 

otherwise 

[23] 

H2  diffusivity  in  gas,  Dh9  (m2  s-1) 

l-lxlO-4(3&)3/20 

0 

[24] 

O2  diffusivity  in  gas,  Do9  (m2  s-1) 

3-2  x  Hr3  (£)3/2(j 

9 

[24] 

H2  diffusivity  in  membrane,  D (m2  s-1) 

2.59  x  10“10 

[3] 

O2  diffusivity  in  membrane,  (m2  s-1) 

1.22  x  1CT10 

[3] 

H2O  diffusivity  in  gas,  Dh9o  (m2  s_1) 

7.35  x  10-5  (  3j3  )  '  ( 

i) 

[24] 

Table  4 


Basic  operating  conditions 


Description 

Symbol 

Value 

Cell  temperature 

Tccll 

80  °C 

Operating  cell  voltage 

Cecil 

— 

Pressure  at  gas  channel  inlet  of  anode  side 

Pa 

2.0  atm 

Hydrogen  mole  fraction  at  anode  gas  channel 

CH2,in 

— 

inlet 

Water  vapour  mole  fraction  at  anode  gas 

CW,in 

— 

channel  inlet 

Stoichiometry  ratio  in  anode  gas  channel 

2.0 

Pressure  at  gas  channel  inlet  of  cathode  side 

Pc 

2.0  atm 

Oxygen  mole  fraction  at  cathode  gas  channel 

c ^ 

C02,in 

— 

inlet 

Nitrogen  mole  fraction  at  cathode  gas 

rc 

CN2,m 

— 

channel  inlet 

Water  vapour  mole  fraction  at  cathode  gas 

CW,in 

— 

channel  inlet 

Stoichiometry  ratio  in  cathode  gas  channel 

2.0 

and  Srinivasan  [21]  that  was  intended  to  apply  an  internally 
humidified  mechanism  to  a  portable  power  application  by 
controlling  the  water  balance.  The  calculated  curves  show 
good  agreement  with  the  measured  experimental  data  for  the 
polarization  curve,  and  also  predict  the  general  trends  of  the 

Table  5 


Design  parameters 


Description 

Symbol 

Value 

Cell/electrode  length 

L 

7.36  cm 

Gas  channel  height 

H 

0.159  cm 

Gas  channel  width 

W 

0.254  cm 

Current  collector  width 

S 

0.082  cm 

Anode  GDL  thickness 

GDLa 

0.02  cm 

Porosity  of  anode  GDL 

£a 

0.4 

Anode  catalyst  layer  thickness 

CLa 

0.001  cm 

Porosity  of  anode  catalyst  layer 

£CLa 

0.112 

Membrane  thickness 

XMW 

Cathode  catalyst  layer  thickness 

CLc 

0.001  cm 

Porosity  of  cathode  catalyst  layer 

£CLc 

0.112 

Porosity  of  cathode  GDL 

£c 

0.4 

Cathode  GDL  thickness 

GDLC 

0.02  cm 

removal  rate  of  water  collected  at  the  anode  outlet  as  a  func¬ 
tion  of  the  hydrogen  flow  stoichiometry. 

Having  validated  the  present  CFCD  model  in  cases  with¬ 
out  external  humidification,  the  following  discussion  involves 
the  general  features  of  water  transport  and  parametric  studies 
to  gain  a  better  understanding  of  the  optimum  water  bal¬ 
ance  over  a  fuel  cell  by  varying  such  parameters  as  the  anode 


1  2  3  4  5  6  7 

Hydrogen  Stoichiometry 
(b) 

Fig.  5.  (a)  Experimental  and  computed  polarization  curves  for  fuel  cell  with¬ 
out  external  humidification  and  (b)  percent  of  water  removed  through  anode 
channel  outlet  as  function  of  hydrogen  stoichiometry  for  problem  described 
in  Fig.  7,  and  comparison  of  experimental  with  numerical  results. 
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(a)  Co-flow 


(b)  Counter-flow 


Fig.  6.  Water  molar  concentration  distribution  in  fuel  cell:  ALCL,  Nation® 


1 1 1,  0.6  V,  and  z  =  Ht. 2.  Left  channel  is  for  anode  and  right  channel  for  cathode. 


and  cathode  flow  directions,  membrane  thickness,  and  inlet 
humidity  values. 

4.2.  Anode  and  cathode  flow  arrangement 

The  predicted  two-dimensional  water  concentration  con¬ 
tours  on  the  x-y  plane  at  z  —  HU  are  shown  in  Fig.  6  for  the 
ALCL  (RH  =  20%)  case  in  the  co-flow  and  the  counter-flow 
designs  of  the  anode  and  cathode,  respectively.  In  the  co-flow 
channel,  water  is  transported  from  the  cathode  to  the  anode 
along  the  entire  gas  channel  (i.e.,  one-way  water  flux),  while 
in  the  counter-flow  design,  water  crosses  through  the  mem¬ 
brane  from  anode  to  cathode  up  to  the  mid-length  along  the 
channel  (i.e.  y  =  U 2),  and  then  water  goes  back  from  cathode 
to  anode  above  the  mid-length  (i.e.,  two-way  water  flux)  as 
shown  in  Fig.  6. 

The  internal  circulation  of  water  resulting  from  the  two- 
way  water  flux  through  the  membrane  is  a  unique  feature 
of  the  counter-flow  design  and  causes  the  electrolyte  mem¬ 
brane  hold  more  water  for  better  ionic  conductivity.  This  is 
evident  in  Fig.  7,  which  shows  the  calculated  water  content 
profiles  across  the  membrane  thickness  at  three  representative 
y-locations  along  the  gas  channel.  In  the  co-flow  design,  the 
water  content  increases  monotonously  with  the  axial  direc¬ 
tion.  Due,  however,  to  the  internal  circulation  of  water  flux 


created  by  the  counter-flow  design,  the  water  content  reaches 
a  maximum  value  in  the  mid-portion  of  the  fuel  cell  (i.e., 
y  =  L/ 2),  while  at  both  ends  of  the  fuel  cell  the  membrane 
contains  less  water,  as  can  be  seen  in  Figs.  6(b)  and  7. 

There  is  an  increase  of  at  least  30%  in  the  predicted  cell 
current  density  for  the  counter-flow  design  at  Vceii  =  0.6  V. 
This  is  ascribed  to  the  internal  circulation  pattern  of  water 
flow  over  the  entire  cell  domain  with  a  thin  membrane  (e.g., 
Nation®  1 1 1).  By  contrast,  there  is  a  negligible  effect  of  the 


Fig.  7.  Water  content  profiles  at  selected  j-points:  Nation®  111,  ALCL, 
0.6  V,  and  z  =  H!2. 
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Table  6 


The  effect  of  operating  parameters  on  the  overall  cell  performance 


Kell  (V) 

(Nafion® 

117)  (/avg,  A  cm  2) 

(Nafion® 

111)  (/avg,  A  cm  2) 

ADCD 

ALCL  (RH  =  20%) 

ADCD 

ALCL  (RH  = 

20%) 

Co-flow 

Counter-flow 

Co-flow  Counter-flow 

Co-flow 

Counter-flow 

Co-flow 

Counter- flow 

0.4 

0.107 

0.108 

0.260  0.271 

1.251 

1.384 

1.337 

1.352 

0.6 

0.037 

0.037 

0.126  0.127 

0.399 

0.590 

0.799 

1.069 

flow  direction  on  the  cell  performance  with  Nation®  117 
when  dry  gases  or  the  low  humidification  of  the  inlet  gases 
are  supplied  as  shown  in  Table  6. 

4.3.  Effect  of  membrane  thickness 

The  predicted  water  concentration  contours  for  Nation® 
117  and  Nation®  112  membranes  on  the  x-y  plane  at  z  =  H/2 
show  very  similar  patterns  of  water  transport  under  the  full 
humidification  of  inlet  gases  in  the  counter-flow  design,  as 
illustrated  in  Fig.  8.  A  greater  production  rate  of  water  is 
expected  on  the  cathode  side  of  Nation®  112  membrane, 
due  to  the  lower  ionic  resistance  of  the  membrane  compared 
with  that  of  the  thicker  Nation®  117  membrane.  Through  the 
anode  gas  channel,  the  enhanced  back-diffusion  of  water  in  a 
thin  membrane  offsets  the  water  loss  by  the  electro-osmotic 


drag,  as  shown  in  Fig.  8.  The  predicted  water  content  pro¬ 
files  across  the  normalized  membrane  thickness  are  given 
in  Fig.  9,  in  which  zero  refers  to  the  anode  vertical  line 
membrane  interface  and  unity  for  the  cathode  vertical  line 
membrane  interface.  The  water  content  level  for  Nation® 
117  and  Nation®  1 12  at  the  cathode  vertical  line  membrane 
interface  is  almost  the  same.  On  the  other  hand,  there  is  a 
large  difference  in  the  water  content  at  the  anode  vertical 
line  membrane  interface,  which  demonstrates  the  effect  of 
the  membrane  thickness  on  the  pack-diffusion  of  water. 

The  predicted  local  current  density  contours  on  the  y-z 
plane  for  Nation®  117  and  Nation®  112  are  presented  in 
Fig.  10  under  the  same  operating  conditions  as  above.  A 
thinner  membrane  yields  more  current  due  to  its  lower  ionic 
resistance.  In  addition,  the  improved  water  back- transport 
for  the  thinner  Nation®  112  membrane  apparently  has  a 


Cj  ,20  [mol/m3] 
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Fig.  8.  Water  molar  concentration  over  fuel  cell  domain:  counter-flow,  AFCF,  0.6  V,  and  z  =  HU. 
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Fig.  9.  Water  content  profiles  at  selected  y-points:  counter- flow,  AFCF,  0.6  V, 
and  z  =  H/2. 

positive  effect  on  the  cell  performance,  as  shown  in  Fig.  10. 
As  discussed  above,  well-hydrated  membrane  shows  good 
ionic  conductivity  which  is  a  function  of  water  content  in 
the  electrolyte  phase.  For  Nation®  117,  the  local  current 
density  has  an  almost  flat  profile  over  the  entire  reaction 
area  and  is  low  compared  with  Nation®  1 12  at  the  same  cell 
voltage,  Veen  =  0.6  V.  In  general,  a  thick  membrane  always 


(b)  Nation®  1 12 


Fig.  10.  Local  current  density  distribution:  counter-flow,  AFCF  at  0.6  V. 


produces  a  more  uniform  current  density  distribution.  By 
contrast,  in  a  thin  membrane,  higher  local  current  densities 
are  predicted  within  the  catalytic  reaction  zone  parallel  to  the 
gas  chamber,  which  signifies  a  large  oxidant  consumption. 
Therefore,  a  lower  available  amount  of  oxidant  results 
in  lower  current  densities  at  the  reaction  area  underneath 
the  current-collector  rib  in  Fig.  10(b)  than  the  minimum 
local  current  density  in  a  thick  membrane,  as  shown  in 
Fig.  10(a). 

4.4.  Effect  of  inlet  gas  humidity 

The  predicted  two-dimensional  contours  of  the  water  con¬ 
centration  on  the  x-y  plane  at  z  =  HI 2  over  the  entire  fuel 
cell  of  Nafion®  1 1 1  at  Vceii  =  0.6  V  for  ADCD  and  ALCL 
(RH  =  20%)  cases  are  presented  in  Fig.  1 1 .  Under  dry  hydro¬ 
gen  and  dry  air  inlet  conditions,  it  can  be  seen  from  Fig.  1 1  (a) 
that  the  water  concentration  remains  below  the  saturation 
level  (i.e.,  16. 14  mol  m-3)  and  results  in  low  ionic  conductiv¬ 
ity  of  the  electrolyte.  Conversely,  adding  a  slight  amount  of 
water  vapour  (i.e.,  RH  =  20%)  helps  the  PEMFC  retain  a  con¬ 
siderable  amount  of  water  over  the  entire  cell  geometry,  as 
illustrated,  in  Fig.  11(b).  The  maximum  water  concentration 
at  the  cathode  vertical  line  membrane  interface  in  Fig.  11(b) 
becomes  greater  than  the  saturation  water  concentration,  i.e., 
16. 14  mol  m-3,  at  80  °C  under  2  atm. 

Computed  profiles  of  water  content  are  presented  in 
Fig.  12  for  selected  y-locations  along  the  gas  channel  under 
the  same  operating  conditions  as  for  Fig.  11.  In  the  ADCD 
case,  the  water  content  A  has  a  constantly  low  value  (i.e., 
approximately  A  =  3)  across  the  membrane  due  to  its  small 
thickness  (Nafion®  111).  The  low  current  output  due  to 
the  large  ionic  resistance  of  the  dry  membrane,  and  there¬ 
fore  an  insufficient  water  production  rate,  cannot  produce  a 
high  water  content  in  the  membrane.  For  ALCL  (RH  =  20%) 
case,  however,  the  water  content  level  on  the  cathode  mem¬ 
brane  interface  exceeds  the  saturation  value  (i.e.,  A  =  14)  and 
there  is  also  a  significant  amount  of  water  at  the  anode  ver¬ 
tical  line  membrane  interface  in  Fig.  12,  which  suggests 
that  water  diffuses  smoothly  back  from  the  cathode  to  the 
anode. 

The  predicted  local  current  density  contours  on  the  y-z 
plane  for  Nafion®  1 1 1  in  the  counter-flow  design  are  shown 
in  Fig.  13  for  operating  conditions  of  Vceii  =  0.6  V,  80  °C,  and 
2  atm  for  ADCD  and  ALCL  (RH  =  20%)  cases,  respectively. 
As  discussed  earlier,  the  thinner  Nafion®  111  membrane  gen¬ 
erally  enhances  the  back-diffusion  of  water.  The  dry  fuel  and 
dry  air  inlet  streams  are  not  enough,  however,  to  cause  a  suf¬ 
ficient  water  concentration  in  the  cell  and  thus  a  conductive 
membrane  for  high  rate  of  water  production.  The  end  result 
is  a  low  cell  current  density.  On  the  contrary,  slight  humidifi¬ 
cation  (e.g.,  RH  =  20%)  of  the  inlet  fuel  and  air  gases  helps  to 
induce  a  positive  chain-effect,  namely,  it  increases  the  water 
concentration  level  so  as  to  provide  good  membrane  con¬ 
ductivity,  improves  the  current  output,  and  intensifies  water 
production  within  the  cell. 
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Fig.  11.  Water  molar  concentration:  Nation®  111,  counter-flow,  0.6  V,  and  z  =  /7/2. 


4.5.  High  performance  regime  of  low -humidity  operation 

Calculated  water  content  profiles  in  the  membrane  are 
given  in  Fig.  14  for  selected  ^-locations  along  the  gas  chan¬ 
nel  of  Nation®  112  at  z  =  HI 2  and  Well  =  0.6  V  for  vari¬ 
ous  humidification  cases,  viz.,  RH  =  20,  50  and  100%.  The 
lowest  humidification  of  the  inlet  gas  streams  (RH  =  20%) 
contains  a  significantly  low  amount  of  water  in  the  mem¬ 
brane,  which  results  in  considerable  ionic  resistance  across 


0.0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1.0 

Normalized  Membrane  Thickness 

Fig.  12.  Water  content  profiles  at  selected  j-points:  counter-flow,  Nation® 
1 1 1,  0.6  V,  and  z  =  HU. 


the  ionomer  membrane.  The  low  cell  current  density,  owing 
to  this  large  ionic  resistance,  cannot  produce  sufficient  water 
by  the  cathode  catalytic  reaction.  With  moderate  humidifi¬ 
cation  (RH  =  50%),  however,  the  water  content  level  in  the 
membrane  is  equivalent  to  that  of  full  humidification  of  the 
inlet  gases  (RH=  100%).  The  water  content  at  the  cathode 
vertical  membrane  interface  reaches  the  saturation  value  (i.e., 
X  =  14)  and  there  is  also  a  significant  amount  of  water  at  the 
anode  vertical  membrane  interface,  as  shown  in  Fig.  14.  This 
results  in  a  high  cell  current  density  comparable  with  that  for 
full  humidification  (RH=  100%). 

4.6.  Polarization  curves 

Computed  performance  curves  for  the  co-  and  counter¬ 
flow  of  Nation®  112  under  the  basic  operating  conditions 
are  presented  in  Fig.  15.  Three  distinctive  characteristics  are 
worth  noting.  First,  with  anode  humidification  only  (AFCD), 
it  is  possible  to  achieve  a  level  of  performance  as  much  as  that 
obtained  with  full  humidification  of  both  anode  and  cathode 
gases  (AFCF).  This  finding  is  of  significant  practical  impor¬ 
tance  because  it  is  easier  to  humidify  the  anode  gas,  due  to 
its  much  smaller  flow  rate,  than  the  cathode  air. 

Second,  either  dry  supply  (ADCD)  or  low  humidification 
(ALCL)  of  inlet  flows  yield  significantly  lower  performance 
than  the  AFCF  and  AFCD  cases  at  high  cell  voltage  of 
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Fig.  13.  Local  current  density  distribution  of  Nafion0  111:  counter-flow  at  ,  ,  .  _ 

_  Fig.  15.  Polarization  and  power  density  curves  for  Nation^  1 12  membrane. 

U.6  V. 


Vceii  =  0. 6V.  As  mentioned  above,  however,  the  performance 
is  much  improved  as  the  cell  voltage  is  reduced.  The  third 
point  is  evident  from  Fig.  15(b),  that  is,  dry  hydrogen  and 
dry  air  feed  streams  (ADCD)  in  the  co-flow  design  cannot 
sustain  a  high  current  density  at  even  low  cell  voltages.  This 
suggests  that  it  would  be  more  suitable  to  adopt  a  counter¬ 
flow  design  to  achieve  the  high  current  density  comparable 


Fig.  14.  Water  content  profiles  at  selected  y-points:  counter-flow,  Nafion® 
112,  0.6  V,  and  z  =  M2. 


with  the  full  humidification  case  at  low  cell  voltage  with  a  loss 
of  less  than  15%  of  the  maximum  power  density,  as  shown 
in  Fig.  15(a). 

5.  Conclusions 

A  unified  water  transport  CFCD  model  based  on  the 
single-domain  approach  has  been  developed  to  account  for 
various  water  transport  modes  in  the  electrolyte  and  to  pro¬ 
vide  a  convenient  tool  to  assess  water  distribution  in  PEMFCs 
that  have  different  humidification  strategies.  The  effects  of 
flow  arrangement,  membrane  thickness,  and  inlet  gas  humid¬ 
ity  on  fuel  cell  performance  have  been  analyzed  to  eluci¬ 
date  the  water  transport  characteristics.  It  is  demonstrated 
that  thicker  membranes  yield  not  only  higher  ohmic  resis¬ 
tances  but  also  less  back-diffusion  flux  of  water  across  the 
membrane.  It  is  also  found  that  low  humidification  of  inlet 
gases  can  produce  a  cell  current  density  equivalent  to  full 
humidification  at  low  cell  voltages  and/or  with  a  counter-flow 
design.  The  latter  design  creates  internal  water  circulation, 
which  facilitates  membrane  hydration.  Finally,  a  fuel  cell  unit 
using  dry  inlet  gases  could  be  an  alternative  power  source 
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for  portable  electronics  provided  that  the  internal  water  cir¬ 
culation  can  be  utilized  to  eliminate  need  for  any  external 
humidification. 
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